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Land snails move via adhesive locomotion. Through muscular contraction and expansion of their 
foot, they transmit waves of shear stress through a thin layer of mucus onto a solid substrate. 
Since a free surface cannot support shear stress, adhesive locomotion is not a viable propulsion 
mechanism for water snails that travel inverted beneath the free surface. Nevertheless, the motion 
of the freshwater snail, Sorbeoconcha physidae, is reminiscent of that of its terrestrial counterparts, 
being generated by the undulation of the snail foot that is separated from the free surface by a thin 
layer of mucus. Here, a lubrication model is used to describe the mucus flow in the limit of small 
amplitude interfacial deformations. By assuming the shape of the snail foot to be a traveling sine 
wave and the mucus to be Newtonian, an evolution equation for the interface shape is obtained 
and the resulting propulsive force on the snail is calculated. This propulsive force is found to 
be non-zero for moderate values of Capillary number but vanishes in the limits of high and low 
Capillary numbers. Physically, this force arises because the snail's foot deforms the free surface, 
thereby generating curvature pressures and lubrication flows inside the mucus layer that couple to 
the topography of the foot. 

I. INTRODUCTION 

Engineers often look to nature's wide variety of locomotion strategies to inspire new inventions and robotic devices 
[H; S; S 0] • More generally, scientists across all disciplines are interested in understanding the physical mechanisms 
behind different styles of biolocomotion. The mechanism of terrestrial snail locomotion has been investigated and 
elucidated over the last couple of decades; conversely, the propulsion of water snails that crawl beneath a free surface 
has yet to be considered. The purpose of this paper is to propose a propulsive mechanism for water snails. 

Gastropod locomotion has been of scientific interest for more than a century @, @, 01 . Three distinct modes of 
locomotion have been examined: ciliary motion, pedal waves, and swimming. Ciliary locomotion, characterized by 
the beating of large arrays of cilia on the animal's foot, is usually distinguished from pedal waves by indirect means, 
such as lack of visible muscle undulation, uniform adherence of the foot to the substrate, and a uniform gliding of the 
snail body This particular type of locomotion is mostly employed by various marine and freshwater snails. 

A significant effort has gone towards understanding pedal wave locomotion by terrestrial snails. Lissmann 0, 
was a pioneer in constructing a mechanistic model of the snail foot undergoing such locomotion. In Ref. Q , he studied 
three species of terrestrial snails {Helix, Haliotis, and Pomatias), all of which use waves of contraction that propagate 
in the direction of their motion ( "direct waves" ) . Waves traveling in the opposite direction ( "retrograde waves" ) were 



examined by Jones and Trueman . A vital insight was later provided by Denny 12[ , who turned the focus of study 
from the snail foot to the properties of the pedal mucus. Pedal mucus has a finite yield stress that allows it to act 
as an adhesive under small strains and to flow like a viscous liquid beyond its yield point. Thus, the snail is able to 
create regions of flow in the mucus by locally shearing it while the rest of the mucus is effectively glued to the solid 
substrate; these regions (or shear waves) propagate along the length of the foot, enabling the snail to move. D enny 



used the nonlinear nature of the mucus to rationalize the locomotion of Ariolimax columbianus, a terrestrial slug 131. 
Recently, it was found that mucus with shear-thinning properties results in energetically favorable locomotion 14| . 
which is well supported by the experimental studies of the mucus properties. A detailed investigation of the rheology 
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FIG. 1: Snail (Sorbeoconcha physidae) crawling smoothly underneath the water surface while the surface deforms. Note the 
surface deflection associated with the undulatory waves propagating from nose to tail along its foot. Photo courtesy of David 
Hu and Brian Chan (MIT). 



of mucus was conducted by Ewoldt et. al [la ], who also tested different synthetic slimes. The versatility of snail 
locomotion has also inspired Chan et. al to build a robotic snail, the first mechanical device to utilize the nonlinear 
properties of this synthetic mucus. 

Terrestrial snails that employ adhesive locomotion are only a fraction of species in the class of gastropods. Tra- 
ditionally, gastropods have been divided into four subclasses: prosobranchia, opisthobranchia, gymnomorpha, and 
pulmonata, the latter of which consists primarily of terrestrial snails. Many species have not been thoroughly inves- 
tigated but exhibit interesting locomotive behavior. For instance, opisthobranchs, that have reduced or absent shells, 

A still more puzzling mode of locomotion is observed in certain species of water snails. 



swim 



16[ or burrow 
In 1910, Brocher ' 



17, 1 



13 ] remarked on water snails that can swim inverted beneath the water surface; since then, other 
qualitative descriptions have been reported. Milnes and Milnes observed the foot of a pond snail "pulsing with 
slow waves of movement from aft to fore along its length," suggesting that direct waves are employed for propulsion. 
The presence of a trail of snail mucus was also reported. Goldacre [2 II ] measured the surface tension of this thin trailing 
film to be approximately 10 dynes/cm; he also remarked that the creature was "grasping the film" as evidenced by 
the film's being pushed sideways as the snail advanced. Deliagina and Orlovsky 22j made similar observations while 
studying feeding patterns of Planorbis corneus. This particular freshwater snail crawls at about 15mm/s, a speed 
comparable to that on land, while the cilia apparent on the organism's sole "beat intensely." Cilia-aided crawling 
beneath a free surface was observed on marine snails as early as 1919. Copcland 23[ concluded that the locomotion 
of Alectrion trivittata, that crawls upside down on the surface, relied solely on the ciliary action. He conducted a 
similar study on Polinices duplicata and Polinices heros, both of which were observed to use both cilia and muscle 
contraction for locomotion on hard surfaces [24J ] . Only ciliary motion was employed by the young Polinices heros 



when crawling inverted beneath the surface. 

It is therefore clear that freshwater and some marine snails have the striking ability to move beneath a surface that 
is unable to sustain shear stresses. In this paper, we attempt a first quantitative rationalization of these observations. 
We use a simplified model based on the lubrication approximation to show that a free-moving organism located 
underneath a free surface can move using traveling- wave-like deformation of its foot. We first present our observations 
of the propulsion of the freshwater snail, Sorbeoconcha physidae in iflTl We introduce our model based on the 
lubrication approximation in £|III[ and present solutions for small-amplitude motion of the foot. The physical picture 
for the generation of propulsive forces is discussed in ^IVl together with the main conclusions and a summary of the 
simplifying assumptions used in our analysis. 



FIG. 2: A trail of mucus behind the snail crawling upside down beneath the free surface. Photo courtesy of David Hu and 
Brian Chan (MIT). 



II. OBSERVATIONS 



Several common freshwater snails, Sorbeoconcha physidae, were collected from Fresh Pond, Massachusetts. About 
1 cm in length, this particular snail can crawl beneath the water surface at speeds as high as 0.2cm/s, comparable to 
its speed on solid substrates, and perform a 180° turn in 3 seconds. It is rendered neutrally buoyant by trapping air 
in its shell. 

The undulation of the snail foot causes surface deformations with a characteristic wavelength of 1 mm and amplitude 
of 0.2 — 0.3 mm (see Figure [l}. This deformation appears to travel in the opposite direction of the snail motion, 
suggesting the generation of retrograde waves, contrary to observations of Milnes and Milnes ■ Another notable 
feature of water snail propulsion is the presence of a trail of mucus (see Figure [2]). For land snails, this mucus layer 
is typically 10 — 20 /im in thickness |15| : as with land snails its rhcological characteristics may also play a significant 
role in underwater locomotion. Since these water snails are also able to crawl on solid substrates, one might venture 
that their mucus properties do not differ too greatly from those of land snails. 



III. MODEL 
A. Assumptions 

The crawling of water snails beneath the free surface has four distinct physical features: a free surface with finite 
surface tension, <r, a layer of (presumably) non-Newtonian mucus, coupled deformations of the foot and the surface, 
and a matching of the flow inside the mucus to that around the snail. To isolate the critical influence of the first 
feature, we consider in this paper a simplified model system characterized by a Newtonian mucus layer and small 
deformations of the foot and the interface, with hopes of providing physical insight into the propulsion mechanism. 



B. General equations 



Choosing a characteristic velocity U ~ lem/s, a mucus thickness H ~ 20 /im, and a (post-yield) mucus viscosity 
v ru 10~ 2 m 2 /s 15| suggests a Reynolds number of the flow within the mucus layer to be Re = UH/v ~ 10 -5 . Thus, 
we neglect inertia and start with incompressible Stokes equations: 



V-v = 0, v-n = o, 



(la, b) 



4 



where v is the velocity field inside the mucus and II the stress tensor. Normal and tangential stress boundary 
conditions at the surface may be expressed as 

n • II • n = an, t • II • n = 0, (2a, b) 

where n and t denote, respectively, unit vectors normal (outward) and tangent to the free surface, and k — — V • n, 
denotes the curvature of free surface. We limit our attention to the two-dimensional case, for which the interface 
shape is given by y = h(x) (see Figure [3]) and n may be expressed as h^$/(l + h~) 3 / 2 , where the "hat" notation 
denotes dimensional variables. The mucus is assumed to be a Newtonian fluid; thus, the stress tensor is given by 

n = -pi + 2ue, (3) 

where e = |{(Vv) + (Vv) T } is the rate-of-strain tensor. 

In the frame moving with the snail, we assume that the gastropod foot undergoes periodic deformations in the form 
of a traveling wave moving at a speed V w . Our approach is as follows: given the shape of the foot, we solve for the 
shape of the liquid-air interface together with the velocity field in the mucus layer, and then calculate the resulting 
propulsive force on the snail. 



C. Lubrication Analysis 

To eliminate temporal variation, we consider the frame of reference moving at the wave speed V w relative to the 
snail (Fig. [3]). We define V w and the (unknown) snail speed, V s , to be positive when the snail moves in the positive 
x-direction while the wave travels in the opposite direction. 

The thickness of the mucus (tens of microns) is observed to be small relative to a typical wavelength of foot 
deformation (millimeters). Thus, we apply the lubrication approximation 25[ 26, 27, 28| and reduce the governing 



equations based on_£T/A = a<Cl, where H is the characteristic thickness of the mucus film, and A the wavelength 
divided by 2n. More specifically, terms of order a or higher arc discarded in the equations of motion. The governing 
equations and boundary conditions (Eq. Q] - Eq. O are non-dimensionalized using the following set of characteristic 
scales: 

(4a) 
(4b) 
(4c) 
(4d) 

(4e) 



X = 


Ax, 


y = 


Hy, 


(u, v) = 


V w (u,av), 


Vs = 


v w v s , 




fj,XV w 


p = 


— — — p. 


H 2 F 



Based on standard lubrication theory, the equations of motion are reduced to the following 

du dv 
dx dy 

+ - 

dx dy 2 ' 
dy 1 



_ du ^ dv 
n x dy ' 

o =-g+£, <5b) 

= (5c) 



while the boundary conditions may be expressed as 

u = 1, at y = h%, (6a) 

v = 0, at y = hi, (6b) 

du . ,„ . 

0=^-, at y = h 2 , (6c) 
dy 

a 3 

■^rh 2 ,xx = ~P, at y = h 2 , (6d) 
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wave frame 
(b) 

FIG. 3: A close-up view of the mucus and the snail foot undergoing a simple sinusoidal deformation of wavelength 2n\. The 
prescribed shape of the snail foot is denoted as hi; the resultant shape of the free surface, /12, is to be solved for. The known 
constant speed of the wave, V w , is set relative to the snail that is translating with an unknown speed, V s . In laboratory frame 
(a), the wave is moving in the negative Si-direction with V w — V s while the snail is moving in the positive s-direction with V s . 
In the frame moving with the wave (b), the snail body appears to move in the positive S-direction with V w . 



where x and y are horizontal and vertical coordinates of the system while z points out of the page, and u and v are 
velocity components in x- and y-directions, respectively. The shape of the snail foot is prescribed by h\(x) while 
h 2 {x) denotes the unknown free surface shape. Note that Ca = [iV w ja is not a Capillary number in the traditional 
sense since V w is not necessarily the characteristic speed of the flow in the lubrication layer. In order to allow surface 
tension effects to remain relevant in the current problem, the curvature term in Eq. (|6d[) is retained despite being 
multiplied by a 3 ; this is a standard practice in thin film problems with surface tension (see, e.g., Goodwin and Homsy 

0)- 

For convenience, we define a modified Capillary number, Ca = /jV w /a 3 a = Ca/a 3 , so that the normal stress 
condition becomes 

-^h 2 ,xx = ~P, at y = h 2 . (7) 
Ca 

Then by integrating Eq. (|5bp twice with respect to y, and applying necessary boundary conditions, we obtain an 
expression for the velocity field in the mucus layer, 

u(x, y) = JLh 2 , xxx (h 2 y - \y 2 + \h\ - h 2 hA + 1. (8) 
Ca V z 1 J 



The resulting volume flux through the layer, 



Q — h-2 — h\ + -=rh 2x xx 

Ca 

is constant since the mucus thickness does not vary with time in this moving reference frame. 

In order to obtain the motion of the snail, it is necessary to consider the forces acting on the organism. Since its 
motion occurs at low Reynolds numbers, the snail is force-free; hence, the forces from the (internal) mucus flow and 
those from the external flow around the body must sum to zero: 

Pint + Fext = 0. (10) 

More specifically, F ext is equal to — -Fdragex, where -Fdrag is the magnitude of the drag force from the external flow; 
F; nt is the traction caused by the flow in the mucus on the foot of the snail and can be expressed as the integral of 
II • rif, where n/ is outward normal to the foot of the snail. In the lubrication limit, Eq. (fT0|) reduces to 

which is a scalar equation representing force-balance in the x-direction. Here n is the number of waves generated by the 
foot and w is the width of the foot in the ^-direction. Physically, the left hand side of Eq. (fTTj) is the propulsive force 
that arises from the internal flow of the mucus and balances the drag from the external flow. Note that Eqs. (|10j) - (fTTT) 
implicitly neglect the overlap regions between the internal mucus flow and the external flow around the organism; we 
will derive in mil HI the asymptotic limit in which this is a valid assumption. 



(fta -Ml 0^1 +h 2 )+h 1 (h 2 - hi) -hi 



(9) 



dh\ du 
dx dy 



P 



dx = Fdi 



(11) 



D. Solution for small-amplitude motion 



In order to solve the model problem, we consider the following limit for foot deformations. If AH denotes the 
typical amplitude of the foot deformation, we define e = AH/X, and assume it to be small. Note that e is a parameter 
independent of the geometrical aspect ratio, a, as can be seen by considering the case where e = 0; in this limit, the 
dimensionlcss parameter, e, is zero when the foot surface is flat, while a remains finite. We choose the foot shape as 



hi = e since, 

and solve for the associated layer profile, h 2 , order by order as 

h 2 = 1 + +e 2 4 2) +0(e 3 ). 
The resulting expression for the flux is given by 
Q 



l + e[h ( 2 1] - sins) + e 2 h% ] 



( £ h 2 . XX x + £2 h 2 .xxx) 



= 1+ 



(eh 



Ca 

(!) +,2,(2) \ 
2,xxx t t ""2 t xxxJ 



Ca 



i(i + .(tf> 

(l + £ (hf> 



sin ir + eh 



21.(2) 



esmx 



21,(2) 



e z h 



l+e\h¥ 



h^ — — since 



e z h 



81.(2) 



e z h 



(12) 



(13) 



(14) 



where Q = Q(°) + eQ^ + e 2 Q^ + 0(e 3 ). Collecting terms of the same order, the leading order 0(1) simply states 
that — 1, and then at order 0(e) we obtain 
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FIG. 4: Free body diagram of a perfectly periodic mucus layer over one wavelength between nodes j and j + 1. Pressures at 
these nodes, p(27rj) and p(2n(j + 1)), as well as the heights, h2(2nj) and /i2(2-7r(j + 1)), are equal by the periodic boundary 
conditions. Above the mucus layer is open to atmosphere with p a t m set to zero. 



This third order linear ODE has an analytic solution given by 

<*"' = « (1> + A > h£») + (5^1) ™ (^) 

. / x \ . I \/3x \ C cos x + sin x . . 

+ 4» exp ) sin j + c2 + 1 , (16) 

where C = -i^, and Ai, A2, and A3 are unknown constants. For convenience, is set to zero by arbitrarily setting 
Q = 1. Note that Q corresponds to the rate of mucus production by the snail. 



E. Boundary conditions 

The real challenge lies in identifying the three independent boundary conditions required to solve for A\, A 2 , and 
A3. As a logical starting point, we proceed by applying periodic boundary conditions over each wavelength: 

/J2,x(27rj) = /i 2 ,x(27r(j + l)), (17a) 
h 2 , xx (2nj) = h 2 , xx {2ir{j + 1)), (17b) 
h 2 ,xxx{^j) = h 2 , xxx {2-K(j + 1)), (17c) 

where j ranges from to n — 1. In this limit, A\, A 2l and A3 vanish, and these boundary conditions yield no motion 
of the snail, regardless of the value of Ca. Conducting a force balance on the mucus layer over one wavelength (as 
shown in Fig. [4]) offers a simple explanation for this result: the periodic boundary conditions ensure that the pressure 
forces acting on the side control surfaces of the mucus layer precisely cancel. Since the top control surface of the 
mucus is exposed to ambient air pressure, there can be no net force that acts on the bottom control surface. Thus, 
no equal and opposite force acts on the snail foot ([1] in Fig. SJ) , suggesting that no net propulsive force is generated 
under strictly periodic boundary conditions. 

Instead, the three boundary conditions should be selected by the physical constraints on snail locomotion. Since 
the moving gastropod is force- and torque-free, the first two conditions should naturally be 

Y,F y = 0, (18a) 
Sr = 0, (18b) 

where F y refers to forces in the y-direction while r is a torque in the z-direction. In all generality, the sum of all 
forces and torques acting on the snail must vanish at low Reynolds number; the forces and torques from the thin 
film of mucus must therefore balance the forces and torques generated by the external flow around the body and 
those arising from gravity. For this analysis, we assume the snail does not rotate, and that its shape is sufficiently 
symmetric that the external viscous torques and y-force vanish. We also assume the organism to be neutrally-buoyant 
and homogeneous, so the forces and torques due to gravity are zero. As a result, Eqs. (|18ap and (|18b[) only require 
the y-ioTCC and z-torque arising from the thin film to vanish. 
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The final boundary condition arises from consideration of the matching between the internal and the external flow 
around the organism. By symmetry, we expect the swimming speed of the snail to be of order V s ~ e 2 . Since we 
are in the Stokes regime, pressure differences across the moving gastropod should scale linearly with the free stream 
velocity, and occur at order e 2 as well. Thus, the pressure difference between the front and back of the snail is zero 
at 0(e), the order of our formulation, and the third boundary condition becomes: 



p(0) = p(2mr). 



These three boundary conditions allow one to obtain a complete expression for h 2 
dimensionless propulsive force, Fprop at 0(e 2 ): 



(i) 



(19) 

and subsequently solve for the 



F, 



prop 



= 3C 



dhi du 
dx~ + ~fy lv= 



P- 



dx 



2mt 



-h 



(i) 



cosx 



+ h 2xxx( h 2 ~ sinx ) + h 2.Lx 



dx. 



(20) 



(2) 

Note that the integral of h 2 xxx vanishes by the matching pressure boundary condition, which is equivalent to h 2 .xx(Q) = 
h 2 ,xx{^ n7T )- Referring back to the constants in Eq. (|16p . A\, in particular, is a non-trival function of Ca. Hence, unlike 
the strictly periodic boundary condition case, the expression for now contains a non-periodic function that gives 
rise to a non-zero propulsive force. 



F. Crawling speed 

To balance the thin-film propulsive force, it is necessary to evaluate the external drag -Fdrag caused by the motion 
of the snail. For simplicity, the snail is modeled as approximately spherical, with radius R. Although approximate, 
this model yields an order of magnitude approximation for the speed of our model snail. Non-dimensionalizing -Fdrag 
by fJ,V w w, the right hand side of (fTTj) becomes Fdrag ~ 67r/y s //i*, where V s is the snail speed scaled by V w and 
\x* = fx/ Hwateri the viscosity ratio of mucus to water. A correction factor / accounts for the aspherical shape of the 
snail as well as the influence of the free surface on the drag coefficient; for the present analysis, / will be treated as 
known for a given crawler. By combining this scaling with the 0(e 2 ) term in the force on the foot generated by the 
mucus layer, the following expression for V s is obtained: 

2 * 

V 3 « ^Fpro P (Ca,n), (21) 

where -Fp rop , the total propulsive force function, is plotted in Fig. EJa) for different values of n. The exact formula 
for -Fprop is not reproduced in this paper for it is long and not informative for the purpose of this analysis, but is 
straightforward to calculate with symbolic packages. Note that the width of the snail foot, w, is taken to be on the 
same order as R; thus, they drop out of Eq. (|2"T|) . 



G. Results 

Fig. [5] shows that the propulsive force vanishes in the limits of both large and small surface tension. In the limit 
of infinitely large surface tension (Ca — > 0), the interface between air and mucus is undcformable and so is analogous 
to a flat surface that cannot sustain shear stress. A snail would simply slip on such surfaces. The detailed behavior 

. . 3 

for small values of Ca is shown in Fig. OJb). The dimensionless force follows the power-law decay Fp rop ~ Ca for 
decreasing Ca for all values of n. 

With zero surface tension (Ca — * oo), a pressure difference across the interface cannot be sustained and so cannot 
drive the flow within the mucus; hence no propulsive force can be generated in this low surface tension limit either. 

In this case, the force follows for all n the power-law decay |-Fp rop | ~ Ca as shown in Fig. [5jc). 

Note that the propulsive force goes from positive to negative at a moderate value of Ca that, in the case of n = 10, 
is around 0.3. Physically, this implies that the snail switches from retrograde waves to direct waves at this critical Ca. 
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FIG. 5: (a) Dimensionless propulsive force, F plop , normalized by the number of wavelengths, n, as a function of the modified 
Capillary number, Ca = /iV w /a 3 a, where the values of n range from 5 to 30 in the increment of 5. In (b) and (c), the absolute 
value of the dimensionless force, -Fprop, is plotted on a logarithmic scale to show the power-law decay in the limits of Ca — » 
and Ca — > oo respectively. The propulsive force exhibits a Ca decay for large Ca, while it decays as Ca for small Ca. 




icr 2 icr 1 10° io 1 

Ca 



FIG. 6: Absolute magnitudes of components of dimensionless propulsive force due to pressure (solid line) and due to shear 
(dashed line) as a function of Ca for n = 10. (note that the shear force is negative). Hence, the total propulsive force which is 
the sum of these two forces is non-zero only when there is a difference between the two. 



In addition, Fig. |5fa) shows that the propulsive force exhibits two distinct maxima for retrograde and direct waves, 
at values of Ca corresponding to 0.15 and 0.8, respectively. Since the maximum propulsive force for the direct waves 
is higher than that for the retrograde, the direct waves may be a faster mode of locomotion for water snails. This 
points to a possible biological advantage of direct over retrograde waves. 

Fig. [5] quantifies the components of propulsive force due to pressure and shear. It is important to note that the 
force due to shear (dashed line) is negative. In the low Ca limit, these two components precisely cancel, leading to no 
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fcl 



26 28 30 32 34 36 
x 
(a) 




hi 

26 28 30 32 34 36 
X 
(b) 

FIG. 7: Dimensionless pressure (a, dashed line) and shear stress (b, dashed line) within the mucus over two wavelengths for 
n = 10. The single dotted line in both (a) and (b) is the shape of the foot, h± while the solid lines describe the shape of the 
interface, /12 for different values of Ca. Black arrows indicate the direction of increasing Ca. 

motion. When Ca is high, they both vanish. For intermediate values of Ca, a difference in the magnitudes of these 
forces results in a net propulsive force. Note the existence of a finite value of Ca for which the propulsive force reaches 
zero, which is a surprising result of our model. 

Going back to the dynamic boundary condition Eq. ((7|) , one can calculate the pressure and shear stress distribution 
to 0(e) inside the mucus layer. These are plotted in Figure [7] for n = 10 along with the shape of the interface, h% ■ 
In the large Ca limit, the interface shape exactly conforms to the shape of the foot hi; in the small Ca limit, the 
interface becomes flat. At intermediate Ca, there exists an asymmetry in the interface shape, associated with the 
exponential term in Eq. (|16[) . that gives rise to the non-trivial propulsive force. The interface shape, pressure and 
shear stresses are plotted for the first wavelength and the last (corresponding to the front and end of the snail) for 
n = 10 in Fig. [8] As shown in this figure, the mucus thickness at the ends deviate substantially from extrapolated 
periodic values creating an asymmetry between the head and the tail. Thus, although the pressures at the ends of 
the crawler are equal at 0(e), there exists a net 0(e 2 ) pressure force acting on the side control surfaces, owing to 
the a 0(e) difference in thickness of the mucus layer (see Fig. To balance this net force, there has to be a force 
acting on the bottom surface of the mucus layer; thus, there exists an equal and opposite force acting on the foot of 
the snail, corresponding to the propulsive force. 

As suggested by Fig. 03 surface tension is the essential ingredient in this mode of locomotion, and the propulsive 
force vanishes in both limits of asymptotically small (large Ca) and large (small Ca) surface tension. The snail would 
therefore have to tunc the way it deforms its foot to exploit the property of the fluid-air interface. As the foot is 
deformed, it forces a lubrication flow in the mucus layer above and leads to the deformation of the free surface. The 
resulting topography of the free surface, constrained by surface tension, is then exploited by the organism to generate 
a propulsive force. In essence, capillary forces transform the flat free surface into a wavy one capable of sustaining 
tangential stresses and so enabling propulsion. 
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FIG. 8: Dimensionless pressure (dashed line) and the interface shape (solid line) in the front (a) and end (b) of the snail for 
n = 10. The dotted line is the shape of the foot, and black arrows are in the direction of decreasing surface tension. 




FIG. 9: Free body diagram of an asymmetric mucus layer across the foot of the snail. (For simplicity, n = 1 in this diagram.) 
Pressures at the ends, p(0) and p(2nn) are equal by the boundary condition; however, they act over two different mucus 
thicknesses, resulting in a net pressure force. 



H. Matching internal and external flows 

In our analysis, we have neglected the fluid forces on the organism, F matc h, arising from the intermediate matching 
region between the internal (mucus) and the external flows. Since the propulsive force of the snail mostly arises from 
the asymmetric shape of the free surface at the head and tail of the foot, this requires further comment. 

Physically, since we are calculating the internal and external flows separately, both need to be considered. In order 
to first estimate the magnitude of F matc h due to the external flow, we refer to the work by Berdan and Leal who 
studied the motion of a sphere near a deformablc fluid-fluid interface. As an extension of previous work in which the 



interface is assumed to be flat 31, 331, the work in Ref. [3fJ considers the limit of small intcrfacial deformation and 
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its effects on the translating body. Unlike our current analysis, the velocity of the sphere is not governed by the shape 
of the free surface but is fixed as U. The small parameter in this paper, e, reduces to a Capillary number, e = pU /a, 
when gravitational effects are not included. Berdan and Leal showed that in the case of a sphere moving parallel to 
the free surface, the deformation of the interface only has a vertical force contribution at 0(e). In the current analysis, 
Fi nt and F oxt , the forces considered in Eq. (fTU)) . are of 0(e 2 ). Therefore, in order to neglect F mat ch consistently, the 
following condition has to be satisfied: 

e 2 < e 2 , (22) 

which requires one to have a look at how e and e are defined. Since U is V s in our problem, we have e ~ pV s /a. 
Recalling from i jlll CI the Capillary number, Ca, is defined in terms of the wave velocity, V w . Because it has been 
shown that V s = V s /V w scales as e 2 , e can be expressed as 

8^.^^ (23) 
a a 

When one replaces Ca with a 3 Ca and rearranges the terms, the criterion to neglect F matc h in Eq. (j22|) reduces to 

eVCa* < 1. (24) 

Since a and e are both small parameters asymptotically approaching zero in the lubrication analysis in the limit 
of small deformation amplitude, e 2 a 6 <C 1, Eq. (|24|) represents a weak constraint on the validity of neglecting forces 
from the intermediate region. 

The second matching force to consider is that induced by the internal flow. Since there is, in general, a height 
difference between the mucus at the front and the back of the snail, the fluid surface will be distorted at either end to 
match with the flat surface far away. Our work will therefore be valid in the limit where the capillary forces resulting 
from these distortions can be neglected, corresponding to an asymptotic limit which we now characterize. 

The two relevant length scales to consider for matching the distorted fluid interface to the flat free-surface in the 
far-field arc the capillary length, l c ~ y 'a y 1 pg (p is the fluid density), and the width of the snail, it;. The fluid interface 
will be distorted over a length, t, into the fluid, where I ~ min(£ c ,ui). The typical curvature pressure arising from 
surface distortion will be on the order of ~ aSh/i 2 , acting on typical height difference Sh between the free surface 
near the snail and the far-field height of the free surface, and therefore contributes to a force on the snail (per unit 
width) on the order of ~ <r(Sh) 2 /£ 2 . Since Sh ~ eH ~ ea\\ the capillary force is on the order of ~ ae 2 a 2 \ / '£ 2 . 
This force has to be compared with that arising from the external flow, given by p e xtR (again, per unit width) where 
p ex t is the typical magnitude of the pressure outside the organism as it is crawling. Since p ex t ~ pV s /R we have 
PextR ~ P-Vs ~ £ 2 py w . The matching condition becomes therefore ae 2 a 2 A / £ 2 -C e 2 p,V w , which is equivalent to 

R 2 /£ 2 < a Can 2 , (25) 

where we have used the estimate R ~ nX. The second matching condition, Eq. (|25[) . requires that the number of 
wavelengths along the snail's foot, n, be sufficiently large. 

IV. DISCUSSION 

In this paper, we have presented a simplified model of water snail locomotion. The physical picture that emerges is 
the following: the undulation of the snail foot causes normal stresses that deform the interface and drive a lubrication 
flow. The resulting stress distribution couples to the topography of the snail foot, leading to a propulsive force. This 
force vanishes in the limit of Ca — » 0, where the interface is flat, and of Ca — > oo, where the topographies of the 
interface and the snail foot precisely match. A finite propulsive force is obtained for intermediate values of Ca. This 
interplay between the free surface and the snail foot distinguishes water snail locomotion from that of their terrestrial 
counterparts. For the latter, the solid substrate on which the snail crawls is fixed; hence, the shape of the snail foot 
alone determines the pressure and shear stresses generated within the mucus layer. For water snails, however, the 
interface is deformed due to the flow created in the mucus by the foot undulation; the interface, in turn, affects the 
dynamics within the mucus layer, creating pressure and shear stresses that act on the foot. This nonlinear coupling 
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between the foot geometry, surface tension, and dynamics within the mucus layer makes the water snail locomotion 
a less straightforward mode of locomotion. 

A direct analogy exists between the thin film comprising the mucus layer of water snails and those arising in coating 
flows; for example, those used in photolithographic processes to fabricate various electronic c omp onents. This class 
of fluids problem has been well studied, both experimentally 33,|3J,|35| and theoretically 
an example of which includes spin coating. Kalliadasis ct al. 



4l| used lubrication theory to show that in the limit of 



small Ca, the interfacial features become less steep, an effect also captured by our model. Mazouchi and Homsy 42 1 
demonstrated that, in the case of large Capillary number, the shape of the free surface nearly follows the topography. 
In the context of water snail locomotion, we saw in mil Dl that the free surface conforms to the shape of the foot in 
the same limit. 

Our study is only the first step towards a quantitative understanding of gastropod crawling beneath free surfaces. 
It is significant in that we have demonstrated the plausibility of locomotion with the minimal ingredients: Newtonian 
fluids and small amplitude deformations. Nevertheless, outstanding issues remain. In the case of adhesive snail 
locomotion on land, the non-Newtonian properties of snail mucus, such as a finite yield stress and finite elasticity, 
play an essential role [l2T ]. Non-Newtonian mucus is likewise expected to have a significant effect for water snails. 
Furthermore, there need to be more systematic observational studies to identify which water snail species exhibit 
which modes of "inverted crawling" . As reported by Copeland 23|, 24 1 and Deliagina and Orlovsky , some species 
of water snails rely entirely on cilia for propulsion beneath the free surface. If such ciliary motion results in no free 
surface deformation, the physical mechanism examined in this paper is of little relevance, and a closer look at the 
cilia-induced flow is suggested. Alternatively, the non-Newtonian properties of the mucus may prove to be significant 
in this case. Categorizing different species according to their propulsion mechanism of crawling (i.e. cilia versus 
muscle contraction) and the constitutive properties of their mucus would provide a more complete physical picture of 
this intriguing form of locomotion. 
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